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Abstract — Apollo 12 and 15 lunar surface magnetometer data with simultaneous 
lunar orbiting Explorer 35 data are used to plot hysteresis curves for the 

+0 . 02U 

whole moon. From these curves a whole-moon permeability ^ = 1.029 o'oi9 

is calculated. This result in^jlies that the moon is not composed entirely 
of paramagnetic material, but that ferromagnetic material such as free iron 
exists in sufficient amounts to dominate the bulk lunar susceptibility. From 
the magnetic data the ferromagnetic free iron abundance is calculated. Then 
for assumed compositional models of the moon the additional paramagnetic iron 
is determined, yielding total lunar iron content. The calculated abundances 
are as follows: ferromagnetic free iron, 5 - ^ wt. total iron in the 

moon, ^ h \rt io. 
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INTRODUCTION 


In this paper we calculate lunar iron abundance from simultaneous 
magnetic field measurements made by instruments on the lunar surface and 
in orbit near the moon. The total iron abundance is the sum of metallic 
free iron and chemically combined iron, calculated for suitable compositional 
models of the lunar interior. 

Previous estimates of whole-moon iron abundance have generally been 
~ 10^ by weight, often based on meteoritic compositional models, using the 
lunar density of 3-3^ s/cm'^ as a constraint. Urey (1962) reported that the 
moon, if composed of chondritic material of proper density, would have ll-l4'5^ 
iron by weight, depending upon the high pressure phases present in the moon. 
Reynolds and Summers (1969), using a mathematical model involving equations 
of state at high pressure, placed a value of 13fo on lunar iron abundance. 

Later Urey and MacDonald (1971 ) estimated values of 3lo total iron using a model 
in which all the iron is present as FeO, and Q.S 3 I 0 total iron for another model 
in which, the iron is combined in equal amounts of FeO and FeS. Wanke et al.(l973) 
reported a value of 9^, modeling 1-unar composition by that of the Allende chon- 
dritic meteorite. These previous results are generally consistent with the 
whole-moon iron abundance calculated in this paper. 

We calculate lunar iron content from whole-moon magnetic hysteresis curves. 
These curves are plots of total surface magnetic field measured by Apollo 12 
and 15 surface magnetometers simultaneous with the external field measured 
by the lunar orbiting Explorer 35 magnetometer. The slopes of the hysteresis 
curves are used to calculate magnetic permeability of the moon, from which we 
can calculate total iron content in the lunar interior. A paramagnetic mineral 
(olivine or orthopyroxene), combined with ferromagnetic free iron, is used 
with the lunar density constraint to calculate the iron abundance. 
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THEORY 


A lunar hysteresis curve is plotted using data obtained during times 
■when the lunar sphere is magnetized by the earth's magnetic field. The 
total magnetic field B measured at a magnetometer site on the lunar surface 
is, in electromagnetic units, 

B = H + , (1) 

where H is the steady external (terrestrial) magnetizing field and M is the / . 

Csee ^'3' 

magnetization field induced in permeable lunar materiaJ^. Alternately, equation \ 

(l) can be expressed B = |o, H, where the relative permeability |i = 1 + 4nk; 

O 

k is magnetic susceptibility in emu/cm'^. 

The lunar interior is modeled by a homogeneous sphere which has an iron 

composed 

Curie point (T^) isotherm at some depth sphere is assumed to be/ 

of a paramagnetic mineral with free iron distributed throughout. The free iron 
is ferromagnetic in the outer shell where T < T^, and paramagnetic in the core 
where T > T^ . The permeability core therefore corresponds to para- 

magnetic material only, while of the shell combines contributions from 
ferromagnetic iron and a paramgnetic mineral. 

Following the method of Jackson ( 1962 ) the total field on the outer sur- 
face of the sphere can be expressed in vector cou 5 )onent form (see APPENDIX for 
details ) : 

B = H^(1+2F)5 + H„(l-F)y + H_(l-F)z (2) 

where 

( 2 q+l)(n,-l) - \3 (t]-1)(2^ + 1 ) 

F= r (3) 

( 2 q+l)(n, + 2 ) - 2 \^(q-l)(p, - 1 ) 
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Here t) = and are relative permeability of shell and core respec- 

tively (permeability of free space = l); X = R 2 /R^j 'where R^ = radius of the 
core and R^ = radius of the whole sphere. The components of B are expressed in 
the ALSEP coordinate system which has its origin on the lunar surface at a 
magnetometer site. The x-axis is directed radially outward from the lunar 
surface; the y and z axes are tangential to the surface, directed eastward and 
northward, respectively. 

Once field measurements of H and B have been used to determine a value for 
F, iron abundances can be calculated for suitable lunar compositional models 
which must meet the 4^itional constraints imposed on the lunar sphere by the 
overall lunar bulk density (3-3^ g/cm"^) and moment of inertia. The core 
radius R 2 is determined by the iron Curie temperature where, because of internal 
temperature and pressure (Bozorth, 1951), iron changes from^ferromagnetic 
paramagnetic state. It is assxumed that the core material has only paramagnetic 

susceptibility k^ (jig = 1 + ^3tk2) and that the temperature variation of the para- 

•ma^netic 

susceptibility varies according to the Langevin theory: 


k = n m /3K T , 


(4) 


where K is the Boltzmann constant, T is absolute temperature, n is the number 

of 

of ions per gram, and m is the atomic moment, m is ^ the order of a few Bohr 

A 

2 + 

magnetrons e.g., for the Fe ion, m = 5-25 M-g "to 5-53 |J-g (Nagata, I 96 I). 

In the outer shell where temperature T < T , the susceptibility k., = k , + 
"apparent" ^ 

k^^, where k^^ is /ferromagnetic susceptibility and k^^ is paramagnetic suscepti- 
bility. The ferromagnetic component is metallic free iron, assumed to be 
composed of multidomain, noninteracting grains. Furthermore, the expected 
pressures and temperatures in the outer shell are such that the ferromagnetic 
susceptibility of iron will not be substantially altered (Bozorth, 1951; 

r 

Kapitsa, 1955). 
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The measured ferromagnetic susceptibility of the shell material is an 

apparent value which differs from the intrinsic ferromagnetic susceptibility 

of the iron because of self- demagnetization of the iron grains and the volume 

fraction of iron in the shell. The apparent ferromagnetic susceptibility 

k is related to the intrinsic susceptibility k„, according to 
al 11 

P 

\l " 1 + N 

where N is the demagnetization factor and p is the volume fraction of the iron. 
For spherical iron grains W = U Jt/3; experimentally N is found to range between 
3 and U (Nagata, I 96 I). We shall use N = 3-5 in our calculations. 
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EXPERIMENTAL TECHNIQUE 


The basic technique uses simultaneous magnetometer measurements of 
the fields B and H, made while the moon is in the quiet geomagnetic tail / 

region (see Fig, 2). The total surface field B is measured by lunar V '5' ^ 

surface magnetometers (LSM) at the Apollo 12, 15, and l6 sites. The LSM \ 

instrument properties are described in detail by Dyal et al. (1972). Simul- 
taneous measurements of the geomagnetic tail field H are made by the lunar 

orbiting Ejcplorer 35 magnetometer, which orbits the moon with 0.5 

J, 

periselentand 5 R aposelene at a period of 11.5 hours. The Expprer 35 
moon • J.-U • n ^ 

in this volume 

orbit is displayed in a"conpanion paper /■ (Dyal et al., 1973 )> and detailed 
instrument characteristics are outlined by Sonett et al., ( 1967 ). 

1. Lunar induction in the geomagnetic field 

The geomagnetic tail consists of two regions ("lobes") of differing 
magnetic polarity, associated with the earth's dipole field extended down- 
stream in the solar wind. At about 60 earth radii the lunar orbit inter- 
sects the magnetic tail as indicated by the accented portion of the orbit 
shown in Fig, 2, and the moon is in the geomagnetic tail field about four 
days of each lunation. The magnetic and plasma morphology of the geotail 
are especially favorable for observing the steady state induced magnetization field 
of the moon. The ambient field in the geotail is steadier and more uniform 
than in any other magnetic region of the lunar orbit. The geo- 

magnetic field shields the moon from the solar plasma, minimizing electro- 
magnetic interaction effects other than the magnetization field. 

In general the total magnetic field at the lunar surface is a vector 
sum of the external magnetic field measured by Explorer 35; the steady 
remanent field at the LSM site; the poloidal field resulting from transverse 
electric (TE) induction; the toroidal field resulting from transverse 
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magnetic (TM) induction; the magnetization field induced in the permeable 
material within the moon; the diamagnetic field associated with the plasma- 
moon interface currents; and the field associated with the dynamical inter- 
action between plasma flow and the lunar fields, A careful consideration 
of each of these contributions to the total field at the lunar surface is 


necessary to justify the data selection and analysis best suited to determine 
accurately the magnetization fields of the moon in the geomagnetic environment. 

We will see that some of these fields are important to our analysis while 
contributions from others can be neglected. 

The solution given by equations (2)and ^ 3 ) foi" the magnetization of the 

moon assumes a uniform external magnetizing field. Under this assumption 

the induced field is dipolafjhigher order multipoles contribute when 

inhomogeneities are present in the external field. Uniformity of the ambient 

geotail field cannot be guaranteed with the available data. It is possible 

to discriminate against data from nonuniform tail fields by requiring magnetic 

conditions at both LSM and explorer instruments simultaneously to be steady 

over long time intervals . If the geomagnetic field and moon are in relative 

motion with a speed characterized by V ^ U/T, where L is the distance between 

LSM and Explorer 35 magnetometers and T is the time over which the field is 

steady, then any nonimiformities can be detected as magnetic fluduations at 

either instrument. Various criteria on the peak-to-peak fluctuations at both 

magnetometers have been used to select data. A summary of the data selection 

criteria will be presented ir\ a Uier subsection. 

sag 

Fluct-uations of the geotail field drive time- dependent eddy currents in the 

a 

lunar interior. These in turn result in/poloidal magnetic field (TE induction) 
which is’ probably the largest of the electromagnetic induction fields measured 


at the surface. In the limit of low-frequency or small-amplitude driving field 
BB„/3t -*■ 0 and , , 

—E' poloidal 

fluctuations, /induction from the /mode vanishes. Therefore the restrictions 
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placed on the peak-to-peak variations of surface and external fields, used 
to discriminate against nonuniform geotail fields, will simultaneously 
eliminate data obviously contaminated by eddy current fields. Furthermore, 
the data averaiging procedures used in the analysis tend to filter out the 
higher frequencies which are the major contributors in this mode. 

Relative motion of the moon and the geomagnetic field results in a 
motional electric field E = V*H in the lunar reference frame, where V 
is the relative lunar velocity and H is the external field. Cxurrents 
driven by E^ set irp magnetic fields (TM induction) that are toroidal about 
the electric field. Theo§;ically this magnetic field is everywhere tangent 
to the lunar surface so that the ALSEP x-coraponents contain minimal toroidal 
induced fields. 

^4Ti'M 

If the induced magnetization field ^ is predominately dipolar, the 
dipole axis will be aligned with the external field, H, At the 1-unar surface 
the radial and tangential components of the dipole field are (see Appendix) 

= 2acos0 and = orsinG, where 9 is the angle between the dipole axis 
and the radius vector to the surface site, and a is the dipole moment. We 
see that the maximum radial component is along the dipole axis, which for 
the average geomagnetic tail field is along the earth-sun line. As a result 
the average surface dipole field at both Apollo 12 and Apollo 15 LSM sites 
is predominantly in the radial or ALSEP x direction. 

As plasma particles spiral about magnetic field lines, their motion 
induces an opposing field resulting in plasma diamagnetism. The solar wind 
plasma diamagnetism has been detected by Colburn et al. ( 1967 ) and is seen 
as a magnetic field increase (B^) of about 1.5 gammas in magnitude in the 
plasma void on the antisolar side of the moon. This diamagnetic field can 
be expressed (Boyd and Sanderson, I 969 ): 

®d “ • 
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‘w\acj'»i’»iv*At of 

H is .the external vacu\am field, N. is the plasma ion density, and T and T. 

A. ’ 1 e 1 

are the electron and ion ten^jeratures perpendicular to- the magnetic field 
H. The average tail magnetic field is about 10 gammas while the plasma 
density is at least two orders of magnitude below average solar wind values . 

Little is known of the tail plasma ten^jerature^but it is probably 

~ 10^ °K. An upper limit on in the geotail is therefore ~ 0.03 gamma, 
which is far below the LSM resolution of 0.2 gamma. Therefore plasma dia- 
magnetism will be neglected in this analysis. 

The dynamic interaction between solar wind plasma and the lunar surface 
fields have been found to be ^ 16 y and ^ ly at Apollo 12 and 15 sites 
respectively (Dyal et. al., 1973)- Plasma streaming pressure responsible 
for compression of surface fields is ^/2^N^mv where m is the proton mass and 
V is the plasma bulk speed. We assiime a typical speed in the tail is characterized 
by the average magnetopause motions which are about 70 km/sec (Mihalov et al., 
1970 ). Inside the magnetotail lobes the plasma density is below about 

O.Ol/cm"^. Therefore the plasma- field interaction in the geotail has to be at 

3 

least 10 times weaker than in the free streaming solar wind and can be neglected. 

Since remanent surface field compression is negligible, it will be assumed 
that the steady remanent fields at both LSM sites are constant in the tail, 
the importance of which will soon be apparent . In addition these fields are 
assumed to be not measiirable by Explorer 35 at periselene. The results of the 
Apollo subsatellite magnetometer reported by Sharp et al. (1972) confirm this 
assumption. Therefore Explorer 35 measures only the external geomagnetic 
field. 

2. Data processing criteria 

Magnetization effects can be maximized and competing electromagnetic 
induction minimized by applying data selection criteria. It has been pointed 
out that at the Apollo 12 and 15 sites the magnetization field is largest in 
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the radial or ALSEP x-axis. Furthermore, toroidal mode fields are not 
present in the x-coraponent and plasma interaction effects are minimal along 
the X direction. Therefore, only the ALSEP x- components of the LSM and 
Explorer 35 data have been use^t to calculate, magnetization fields. Several 
averaging intervals (6 minutes, 15 minute^ 30 minute^ 1 ho\ir) were used for 
analysis of Apollo 12 data, each with different peak-to-peak criteria, to 
insure external field uniformity and minimize poloidal field contamination. 

Only 60-minute and 3-hour averages, without peak-to-peak criteria, were 
useful for the Apollo 15 data since the simultaneous Explorer 35 data were 
conta[^minated by a high-frequency satellite spin-modulated tone. A typical 
30-minute swath of Apollo 12 data used cn the analysis is shown in Fig. 3- ' 

The steadiness of thesedata is characteristic of much of the geomagnetic 
tail field. 

3. Lunar Hysteresis curves 

Typical lunar "hysteresis" curves are shown in Fig. 4. Apollo 12 and 
15 averages of x components are plotted on the abscissa with the simultaneous 
Explorer 35 averages on the ordinate. These. data were taken from four luna- 
tions of Apollo 12 and two lunations of Apollo 15 data. Bifurcation in the 
distributions is a result of the data selection criteria which eliminate 
data from near the neutral sheet . The normal "S" shape of the hysteresis 
curve degenerates at these low field values to a straight line (EUwood, 1934) 
intersecting the origin since the x axis intercept (x-component of the remanent 
field at both sites) has been subtracted from the data. 

4. Statistical Analysis 

Four lunations of Apollo 12 data and two of Apollo 15 data have been 
processed for each of the averaging times, employing the data selection criteria 
described previously. This has resulted in a total of six separate hysteresis 
curves, one for each averaging interval of both instruments. To each of these 



data sets a least squares straight line has been fitted. , The average and 
standard deviation values of the slopes are 1.019 - 0.012. Both extrema 
are greater than 1.0, implying that the moon, as a whole, acts as a para- 
magnetic or weakly ferromagnetic sphere. 

Scatter in the Apollo 12-Explorer 35 data points of Fig. 4(a) is primarily 
a result of magnetic inhomogeneities between the moon and Explorer 35) small 
contributions from eddy current fields, and instrumental noise in the Apollo 
and Explorer magnetometers. These error sources may introduce small random 
fluctuations into the data which will not substantially affect the slope or 
intercept of the least-squares line. An error source which could affect the 
least-squares slope is a time -dependent heading drift of the Explorer 35 magne- 
tometer. The drift problem, which began after the Apollo 12 and Explorer 35 
data used in Fig. 4(a) were acquired, is therefore applicable only to analysis 
of simultaneous Explorer 35 and Apollo 15 data. The slope of the Apollo 15- 
Explorer 35 hysteresis curve (1.019 - 0.004) shown in Fig. 4(b) agrees well 
with that obtained from the Apollo 12 -Explorer 35 curve (1.019 * 0.013)) 
however, indicating that the Explorer 35 heading error is minimal for the time 
periods used in Fig. 4(b). Furthermore, the agreement in the analysis of data 
from these two Apollo surface sites, which are separated by 1000 km, is a good 
test of our assumption of lunar azimuthal symmetry. 
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RESULTS 


1. Bulk magnetization properties of the lunar sphere 

It is shown in the Theory section that the slope of the hysteresis curve 

at very low field values is related to the magnetic susceptibility of the 

moon. Fig. 5 relates F of equation (3) to the susceptibilities k, and k^ 

= 1 ^ 

of the lunar shell and core. The curve e = k^/kg /corresponds to the case 

of a magnetically homogeneous sphere of permeability jx = 1 + 4-nkj^, 

from which bulk permeability of the moon is calculated to be jx. = 1.029 ^o'oTq * 

Error limits are one standard deviation. -u.019 

^Since both upper and lower limits of permeability are above 1.0, the moon as 

a whole responds magnetically as a paramagnetic or weakly ferromagnetic sphere. 

The field at the surface is in general a function of both shell and core sus- 
ceptibilities; as the susceptibility ratio e = i^^^reases, the core 

susceptibility kg becomes less important magnetically, and the family of 
curves merge for ® ^ 100. 

A further use of the bulk liinar permeability is to place limits on the 
magnetic dipole moment induced in the moon by the geomagnetic field. The dipole 
moment is HR^F (see Appendix); it is oriented parallel to the external field, 

+ 25 

and for H = 10 gammas, the induced dipole moment is of magnitude 5-0 - 5«2 x 10 

gammas-cm^. Russell et al. (1975) using subsatellite magnetic data, have placed 

25 3 

an upper limit on the permanent lunar dipole of I.7 x 10 gammas-cm . While 
interesting to compare, these two results are for different phenomena. The 
averaging techniques performed over many sub satellite orbits to obtain the 
permanent dipole value may effectively average out the magnetization fields 
which change magnit\ide and polarity as a function of the external field H . 
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2. Lunar composition and iron abundance 


In this section we calculate iron abundance in the moon, using various 
compositional models of the lunar interior. The models are constrained by 
magnetic permeability, temperature, pressure, density, and moment of inertia 
determinations . 

The limar moment of inertia is approximately that of a sphere of homo- 
geneous density (Ringwood and Essene, 1970). We will constrain our lunar 

models to be of -uniform density p = 3-3^ gm/cm^ throughout. 

at 

We have selected "two paramagnetic miners to model the lunar interior: 
oli-vine and orthopyroxene. In the moon's outer shell, free iron is assumed to lee 
interspersed throughout the paramagnetic rock matrix. These minerals, 

•with iron content constrained by the overall density of 3 -3^5 have properties 
which are consistent with measured seismic velocity profiles at mantle 
depths (Toksoz et al. , 1972) and geochemical considerations (Ringwood and 

Essene, 1970; Green et al., 1971; Gast, 1972). Furthermore, pyroxenes and 

olivines have been reported to be major mineral components of lunar surface 

Nagata et al. , 1971; Zussman, 1972; - 

fines and rock samples ^Weeks, 1972), with combined iron present as Fe . The 

ferromagnetic conponent of l\inar samples is primarily metallic iron (Wagata 

et al., 1972; Pearce et al., 1971 ) most of which is thought to be native to 

the moon rather than meteoritic in origin (Strangway et al., 1975). 

Magnetization of the l\anar sphere, as measured by the Apollo-Explorer 55 
magnetometer network, imposes a magnetic constraint on the lunar conposition. 
This constraint is determined by the hysteresis-curve slope, represented by the 
value F = O.OO95 - O.OO60 (see equations (2) and (5)); F is in turn related 
to the internal susceptibility of the moon. 

The Curie point of metallic iron will defines the boundary between a homo- 
geneous shell of susceptibility k^ and a homogeneous core of susceptibility 
k^. The boundary is located using the tenperature profiles proposed for the 
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moon by several authors as shown in Fig. 6, Also plotted is the Curie 
temperature of iron which is a function of pressure and therefore of 
depth (Bozorth, 1951). 

For the calculations that follow, three temperature models are used 

'iOar*«V 

in an effort to match hot, and*cold models. For the hot 

model the Curie point is at A = = 0.9. We assume shell and core 

temperatTires are 600 °C and l400 °C, respectively. For the W(Xwv> 
model the shell temperature is 500 °C and that of the core is I 3 OO °C, 
while the boundary is at A = O. 85 . Temperatures are 3 OO °C and 7 OO °C for 
shell and core of the cold model, which has A = 0.7. The magnetic properties 
of the core are fentirely ^determinedj by the paramagnetism of olivine 
(Nagata et al. , 1957) and orthopyroxene (Akimoto et al., 1958) and the 
Langevin tanperature dependence of these minerals (see equation (4)). Small 
amounts of free iron will not substantially change the total effective suscepti- 
bility of the core since iron is paramagnetic above its Ciirie temperature. 

Crustal susceptibility will be determined by the free and combined 
iron under the assumption that the free iron is in noninteractive, multi- 
domain grains whose shape demagnetizing factor is on the average 3-5 (Nagata, 
1961 ). Normally the ferromagnetism of free iron is dominated by a dependence 
on pressure and temperature. However at very low field strengths, such as 
those of the geomagnetic tail, the susceptibility of iron is not strongly 
dependent on temperature below the Curie point (Bozorth, 1951). Uniaxial 
stress on iron changes its susceptibility (Kern, I 96 I); however, hydrostatic 
stress should not affect the susceptibility (Kapitsa, 1955) unless, at great 
pressures, there is a change in volume (Breiner, I 967 ). Therefore, it will be 
assumed that the susceptibility of free iron in the shell is independent of 
pressure and temperature. 




Using the constraints outlined above, forward calculations have been 
made to determine the iron content in the moon. Olivine (lO mole ^ Fe^SiOj^) 
and orthopyroxene (25 mole FeSiO^) minerals are selected using the bulk lunar 
density constraint. For a given temperature model Langevin’s law yields a 
range of susceptibilities for each core composition (see Fig. 7 ) from which 
the paramagnetic iron content of the core is calculated. 

Core susceptibility is then used, along with the magnetically determined 
range of F, to define a range of susceptibility for the shell. Apparent sus- 
ceptibility of the shell is then a function of the free and combined iron in the 
shell. In all cases it is found that the shell susceptibility cannot be 
accounted for by the paramagnetic susceptibility alone; rather, the 

apparent ferromagnetic susceptibility (k^^. ^ dominant: ^ 100. 

Therefore, using equation ( 5 ) the shell free iron abundance is determined to 
accovint for the range of F values (see Fig. 5)- Finally, the paramagnetic 
rock matrix of the shell is chosen so that, containing the free iron, it matches 
the bulk lunar density. Both the ferromagnetic and paramagnetic iron content 
of the shell are then determined. 


From these calculations it has been found that the paramagnetism of any 
olivine and orthopyroxene cannot account for the measured bulk permeability 
of the moon^ contributions from free iron or other highly permeable material 
are necessary. When metallic iron is included in a relatively thin outer shell 
("hot" model where XSO. 9 ), it is difficult to account for the measured bulk 
permeability of the moon while still maintaining the constraint of homogeneous 
density throughout the moon. 

When the iron Curie point is asstuned to be deeper in the moon, i.e., X=0.85 
("warm" model) or X=0.7 ("cold" model), magnetic permeability and density con- 
straints can be easily satisfied by our conpositional models. Ferromagnetic 
iron content in the shell is the same for both olivine and orthopyroxene. 


-15- 



but is a function of thermal model. The warm model requires 2 to 9 'wt.^ free 
iron in the shell, -whereas the cold model requires 1 to 5 wt.^. The total iron 
abundance, on the other hand, is the same for both thermal models, but different 
for the -two mineral compositional models. For olivine (mole % Fe^SiOj^ 5- 10^) 
total iron is 5 to 6 wt.'^; for orthopyroxene (mole % FeSiO^ < 25^), total 
iron is 12 to 13 'wt.’/o. The iron abundance limits reported here all include the 
one- standard-deviation limits of the measured hysteresis cur-ve slope 1.019^0.012. 
Results are summarized in Fig. 8. 



SUMMARY AND CONCLUSIONS 


1. Apollo 12 and 15 lunar surface magnetometer data^ simultaneous with Explorer 

35 magnetometer dataj are used to plot hysteresis curves for the whole moon; 

+0 C2k 

from these curves a whole-moon permeability p,=1.029 is calculated* Both 

error limits constrain the permeability to be greater than 1»0, implying that 
the moon as a whole is paramagnetic or weakly ferromagnetic. 

2. The permeability calculation places a constraint on the composition of the 
limar interior. When the moon is modeled by homogeneous paramagnetic material 
(olivine or orthopyroxene) of density 3-3^ gm/cm^, it is found that these 
purely paramagnetic models have magnetic permeabilities far lower than the 
range determined from magnetometer measurements. We have found no reasonable 
paramagnetic mineral or combination of minerals, with the correct lunar density, 
that would have permeability high enough to be consistent with our measured 
values. From this we infer that the moon must contain some material in the 
ferromagnetic state, such as free metallic iron, in order to accoimt for the 
measured global permeability. 

3. The paramagnetic-ferromagnetic combination is investigated using a two- layer 

lunar model; the boundary between the layers is determined by the location of 

the iron Curie point. Three thermal models, "hot", "warm", and "coldy are , 

used to establish the Ctirie isotherm location. "Warm" and "cold" models fit 

both the density and permeability constraints whereas the "hot" models do not. 

We conclude that the ferromagnetic free iron abundance of the outer shell is 2 to 

yielding an overall range of 5^4 wt.^. 

9 Vt.^ for "warm" andIto5wt.% for "cold" models The free iron content is 

dependent upon the volume of lunar material below the Curie tenperature and 

independent of the composition of the nonferromagnetic portion of the moon. 

4. For two compositional models the total iron abundance is calculated (see 

lTicMin9 

Fig, 8). ^both the olivine and orthopyroxene models, the total iron is 
independent of the thermal model used and is within the limits of 9±4 wt.^o. 
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APPENDIX 


Consider a radially inhomogeneous^o-iayer permeable sphere in an initially 

uniform magnetic field In the absence ofctirrents H = -V5 and since B = pB, 

at any point V-H = 0. Therefore $ satisfies Laplace's equation consistent with 

the continuity of the normal components of B and tangential components of H at 

the spherical boundaries. In the various regions the potentials must be 

00 

V 

- H^r cos ® + 2. ~ i Pj(cos 0) r > 

° t=o r 

CO 

* = \ I B^< T< 

1=0 r 

00 

Y Pj(cos 0) r < Rg 

1=0 

L 

with R^ and the outer radius and permeability, and Rg and pg the inner 
radius and permeability. Using Jackson's (1962) notation conventions we 
match the boundary conditions 


||(R+) = ||(R-) ||(R+) - H(Rg) 

^(Ri) = M-i ^(Ri) M-i ^(Rp = m-2 |?(RP 


All but the t = 1 coefficients vanish leaving the four simultaneous equations 

« - - Y = 


2a + 

11 
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r|p + Y 
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II 
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P^r|p - 2h^y 

- PgR^S =0 


23 - 



We are interested in the induced dipole moment, which is given 

by 


^ (2Tl+l)(^i^-l) - x3(ti-1)(2p,^+1) 

(2Tv^l)(^i^+2) - 2 X^(t)-D(m.^-1) 


where t) = ^- 2/^2 ^ ~ ^2/^1* resulting dipolar field in an (f, 6, i) 

spherical coordinate system is 

B = (2 cos 0 r + sin 6 0) . 

Rl 

We rotate the coordinate system about the r (or ALSEP x) axis so that 0 and $ 
correspond to the ALSEP y and z axes ; then 

B = F (2H X - H y - H z) 

-o X y z ' 


The total field at the surface is the sum of the \aniform external field H 

and the dipole field B : 

— o 

B = (1+2F) X - Hy (1-F) y + (1-F) z . 
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FIGURE CAPTIONS 


Fig. 1. Induced magnetization of the moon. When a two-layer permeable 
sphere is immersed in a uniform external magnetic field H (in our case 
the moon in the steady geomagnetic tail field), a dipolar magnetization 
field M is induced with its dipole moment aligned along the direction of 
H; the total field near the sphere is B = H + UitM. In the case illustrated 
here, > M-g > = 1 (free space), corresponding qualitatively to the 

case where the outer shell (of radius R^) contains both ferromagnetic and 
paramagnetic material, whereas the core (of radius and tenperature above 
the C\jrie point) contains only paramagnetic material. Measurements of B and H 
allow construction of a classical B-H hysteresis curve for the sphere, from 
which permeability and iron abvindance can be calculated. 

Fig. 2. Magnetic environment of the moon during a conplete lunar orbit, 
with emphasis on the geomagnetic tail region. The plane of the lunar orbit 
very nearly coincides with the ecliptic plane of the earth's orbit. The 
earth's permanent dipole field is^ swept back into a cylindrical region known 
as the geomagnetic tail; at the lunar distance the field magnitude is ~ 10 

-4 

gammas or 10 oersteds. Substructure of the tail consists of two "lobes" j 
the vipper or northward lobe has its magnetic field pointing roughly toward 
the earth, whereas the lower lobe field points away from the earth. The 
moon is immersed about four days of each orbit in .the tail; the moon can pass 
throtigh either or both lobes (accented portion of orbit), depending upon the 
characteristics of the particular orbit. 

Fig. 3- Simultaneous Apollo 12 and Ames Explorer 35 magnetometer data from 
a quiet region of the geomagnetic tail, 1970. Vertical scales for the two 

r 

magnetometers differ due to the existance of a 38 ± 3 gamma remanent magnetic 
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field at the Apollo 12 surface site. Components of the fields are expressed 
in the AISEP coordinate system, which has its origin located on the moon's 
surface at the Apollo 12 site; the x-axis is directed radially outward from 
the lunar surface, whereas y and z eoces are tangential to the surface, 
directed eastward and northward, respectively. 

Fig. h. Hysteresis curves for the moon. For each data point the moon is 
immersed in the external magnetizing geomagnetic tail field H (measured 
by Explorer 35) which is steady over 30-niinute or one-hour time intervals; 
the total magnetic induction is B = pH, measured by an Apollo surface mag- 
netometer. Here only the radial (x) conponents are plotted. In this low- 
applied- field regime (~ 10 gamma or lO”^ Oe), the hysteresis curves are 
linear and are fitted by least-squares straight lines, (a) Simultaneous 
Apollo 12 and Explorer data; each data point represents a 30-niinute average 
during quiet times in the geomagnetic tail for fotB?* lunations, 1969-1970. 

(b) Simultaneous 1-hour averages of Apollo 15 and Explorer data for two 
lunj^tions, 1971 - 

Fig. 5. The function F versus susceptibility of the lunar outer shell 
(refer to equations (2) and (3)). This figure is for the case where the 
iron Curie temperature isotherm is at core radius X = Rg/^l “ 0.85, as 
inferred from tenperature models illustrated in Fig. 6. Superimposed on 
the vertical scale are the best value and limits of F calculated from lunar 
hysteresis curve slopes; the calculated range of F determines possible 
ranges for susceptibilities of core and shell for this model. 
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Fig. 6. Ten 5 >erature profiles for the lunar interior published hy various 
authors: 1. McConnell et al. (1967 )j 2, Fricker et al. (196?), 3. Hanks and 

Anderson (1972), k. Dyal and Parkin (1973), and 5- Sonett et al. (1971). 
Superin 5 )osed is the pressure-dependent Curie temperature for iron (Bozorth, 
1951 ) versfls depth in the moon. The Curie point is reached at a position 
R/R ~ 0. &5 for model and ~ 0.7 for a "cold" model. 

Fig. 7" ^fagnetlc susceptibility versus density for two minerals at I 6 OO °K, 
olivine (xFe^SiOi^* (l-x)^^®^^!^) orthopyroxene (xFeSiO^* (l-x)l^SiO^) , 
used to model the lunar interior. For a homogeneous moon composed of either 
of these minerals, the iron silicate concentrations with densities greater 

O 

than p = 3 . 3 ^ gm/cm'^ are forbidden. 

Fig. 8. Iron abundance in the moon. Total iron abundance of the moon is 
calculated for two compositional models of the lunar interior, each with 
two different thermal profiles. Total iron includes both free iron and 
iron combined in the paramagnetic state. The portion of the total iron 
which contributes most strongly to* the measured permeability, ferromagnetic 
free iron, is listed separately. 
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